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Abstract 

Using a family of modified Weibull distributions, encompassing both sub-exponentials and super-expo- 
nentials, to parameterize the marginal distributions of asset returns and their multivariate generalizations 
with Gaussian copulas, we offer exact formulas for the tails of the distribution P(S) of returns S of a 
portfolio of arbitrary composition of these assets. We find that the tail of P(S) is also asymptotically 
a modified Weibull distribution with a characteristic scale \ function of the asset weights with differ- 
ent functional forms depending on the super- or sub-exponential behavior of the marginals and on the 
strength of the dependence between the assets. We then treat in details the problem of risk minimization 
using the Value-at-Risk and Expected-Shortfall which are shown to be (asymptotically) equivalent in this 
framework. 



Introduction 

In recent years, the Value-at-Risk has become one of the most popular risk assessment tool (Duffie and 
Pan 1997, lorion 1997). The infatuation for this particular risk measure probably comes from a variety of 
factors, the most prominent ones being its conceptual simplicity and relevance in addressing the ubiquitous 
large risks often inadequately accounted for by the standard volatility, and from its prominent role in the 
recommendations of the international banking authorities (Basle Commitee on Banking Supervision 1996, 
2001). Moreover, down-side risk measures such as the Value-at-risk seem more in accordance with observed 
behavior of economic agents. For instance, according to prospect theory (Kahneman and Tversky 1979), the 
perception of downward market movements is not the same as upward movements. This may be reflected 
in the so-called leverage effect, first discussed by (Black 1976), who observed that the volatility of a stock 
tends to increase when its price drops (see (Fouque et al. 2000, Campbell, Lo and McKinley 1997, Bekaert 
and Wu 2000, Bouchaud et al. 2001) for reviews and recent works). Thus, it should be more natural to 
consider down-side risk measures like the VaR than the variance traditionally used in portfolio management 
(Markowitz 1959) which does not differentiate between positive and negative change in future wealth. 
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However, the choice of the Value-at-Risk has recently been criticized (Szergo 1999, Danielsson et al. 2001) 
due to its lack of coherence in the sense of Artzner et al. (1999), among other reasons. This deficiency 
leads to several theoretical and practical problems. Indeed, other than the class of elliptical distributions, 
the VaR is not sub-additive (Embrechts et al. 2002a), and may lead to inefficient risk diversification policies 
and to severe problems in the practical implementation of portfolio optimization algorithms (see (Chabaane 
et al. 2002) for a discussion). Alternative have been proposed in terms of Conditional- VaR or Expected- 
Shortfall (Artzner et al. 1999, Acerbi and Tasche 2002, for instance), which enjoy the property of sub- 
additivity. This ensures that they yield coherent portfolio allocations which can be obtained by the simple 
linear optimization algorithm proposed by Rockafellar and Uryasev (2000). 

From a practical standpoint, the estimation of the VaR of a portfolio is a strenuous task, requiring large com- 
putational time leading sometimes to disappointing results lacking accuracy and stability. As a consequence, 
many approximation methods have been proposed (Tasche and Tilibetti 2001, Embrechts et al. 2002b, for 
instance). Empirical models constitute another widely used approach, since they provide a good trade-off 
between speed and accuracy. 

From a general point of view, the parametric determination of the risks and returns associated with a given 
portfolio constituted of N assets is completely embedded in the knowledge of their multivariate distribu- 
tion of returns. Indeed, the dependence between random variables is completely described by their joint 
distribution. This remark entails the two major problems of portfolio theory: 1) the determination of the 
multivariate distribution function of asset returns; 2) the derivation from it of a useful measure of portfolio 
risks, in the goal of analyzing and optimizing portfolios. These objective can be easily reached if one can 
derive an analytical expression of the portfolio returns distribution from the multivariate distribution of asset 
returns. 

In the standard Gaussian framework, the multivariate distribution takes the form of an exponential of minus 
a quadratic form X'Q~ 1 X, where X is the uni-column of asset returns and Q is their covariance matrix. 
The beauty and simplicity of the Gaussian case is that the essentially impossible task of determining a 
large multidimensional function is collapsed onto the very much simpler one of calculating the N(N + 
l)/2 elements of the symmetric covariance matrix. And, by the statibility of the Gaussian distribution, 
the risk is then uniquely and completely embodied by the variance of the portfolio return, which is easily 
determined from the covariance matrix. This is the basis of Markowitz (1959)'s portfolio theory and of the 
CAPM (Sharpe 1964, Lintner 1965, Mossin 1966). The same phenomenon occurs in the stable Paretian 
portfolio analysis derived by (Fama 1965) and generalized to separate positive and negative power law tails 
(Bouchaud et al. 1998). The stability of the distribution of returns is essentiel to bypass the difficult problem 
of determining the decision rules (utility function) of the economic agents since all the risk measures are 
equivalent to a single parameter (the variance in the case of a Gaussian universe). 

However, it is well-known that the empirical distributions of returns are neither Gaussian nor Levy Stable 
(Lux 1996, Gopikrishnan et al. 1998, Gourieroux and Jasiak 1998) and the dependences between assets are 
only imperfectly accounted for by the covariance matrix (Litterman and Winkelmann 1998). It is thus desir- 
able to find alternative parameterizations of multivariate distributions of returns which provide reasonably 
good approximations of the asset returns distribution and which enjoy asymptotic stability properties in the 
tails so as to be relevant for the VaR. 

To this aim, section 1 presents a specific parameterization of the marginal distributions in terms of so-called 
modified Weibull distributions introduced by Sornette et al. (2000b), which are essentially exponential of 
minus a power law. This family of distributions contains both sub-exponential and super-exponentials, 
including the Gaussian law as a special case. It is shown that this parameterization is relevant for modeling 
the distribution of asset returns in both an unconditional and a conditional framework. The dependence 
structure between the asset is described by a Gaussian copula which allows us to describe several degrees of 
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dependence: from independence to comonotonicity. The relevance of the Gaussian copula has been put in 
light by several recent studies (Sornette et al. 2000a, Sornette et al. 2000b, Malevergne and Sornette 2001, 
Malevergne and Sornette 2002c). 

In section 2, we use the multivariate construction based on (i) the modified Weibull marginal distributions 
and (ii) the Gaussian copula to derive the asymptotic analytical form of the tail of the distribution of returns 
of a portfolio composed of an arbitrary combination of these assets. In the case where individual asset 
returns have modified- Weibull distributions, we show that the tail of the distribution of portfolio returns S 
is asymptotically of the same form but with a characteristic scale \ function of the asset weights taking 
different functional forms depending on the super- or sub-exponential behavior of the marginals and on the 
strength of the dependence between the assets. Thus, this particular class of modified-Weibull distributions 
enjoys (asymptotically) the same stability properties as the Gaussian or Levy distributions. The dependence 
properties are shown to be embodied in the N(N + l)/2 elements of a non-linear co variance matrix and the 
individual risk of each assets are quantified by the sub- or super-exponential behavior of the marginals. 

Section 3 then uses this non-Gaussian nonlinear dependence framework to estimate the Value-at-Risk (VaR) 
and the Expected-Shortfall. As in the Gaussian framework, the VaR and the Expected-Shortfall are (asymp- 
totically) controlled only by the non-linear covariance matrix, leading to their equivalence. More generally, 
any risk measure based on the (sufficiently far) tail of the distribution of the portfolio returns are equivalent 
since they can be expressed as a function of the non-linear covariance matrix and the weights of the assets 
only. 

Section 4 uses this set of results to offer an approach to portfolio optimization based on the asymptotic 
form of the tail of the distribution of portfolio returns. When possible, we give the analytical formulas of 
the explicit composition of the optimal portfolio or suggest the use of reliable algorithms when numerical 
calculation is needed. 

Section 5 concludes. 

Before proceeding with the presentation of our results, we set the notations to derive the basic problem 
addressed in this paper, namely to study the distribution of the sum of weighted random variables with given 
marginal distributions and dependence. Consider a portfolio with m shares of asset i of price pi(0) at time 
t = whose initial wealth is 

N 

W(0) = ^n m (0) . (1) 
i=l 

A time r later, the wealth has become W(t) = Yl%=i n iPi( T ) an d the wealth variation is 

N ( \ _ (f)\ N 

S T W = W(t) - W(0) = V m P i(0) Pt{T) , n f W = W(0) V WiXi (t, r), (2) 

where 

niPi(0) 

Wi = (3) 

Ej=i n iPi(°) 

is the fraction in capital invested in the zth asset at time and the return Xi(t, r) between time t — r and t of 
asset i is defined as: 

Xl (t,r) = P ^-^ t - T) . (4) 
Pi{t ~ r) 

Using the definition ©, this justifies us to write the return S T of the portfolio over a time interval r as the 
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weighted sum of the returns rj(r) of the assets i = 1, N over the time interval r 



5 T W 
W(0) 



N 



St 



y^wi Xi(r) . 



(5) 



i=l 



In the sequel, we shall thus consider the asset returns Xj as the fundamental variables and study their 
aggregation properties, namely how the distribution of portfolio return equal to their weighted sum derives 
for their multivariable distribution. We shall consider a single time scale r which can be chosen arbitrarily, 
say equal to one day. We shall thus drop the dependence on r, understanding implicitly that all our results 
hold for returns estimated over time step t. 

1 Definitions and important concepts 
1.1 The modified Weibull distributions 

We will consider a class of distributions with fat tails but decaying faster than any power law. Such possi- 
ble behavior for assets returns distributions have been suggested to be relevant by several empirical works 
(Mantegna and Stanley 1995, Gourieroux and Jasiak 1998, Malevergne et al. 2002) and has also been as- 
serted to provide a convenient and flexible parameterization of many phenomena found in nature and in 
the social sciences (Laherre and Sornette 1998). In all the following, we will use the parameterization 
introduced by Sornette et al. (2000b) and define the modified-Weibull distributions: 

Definition 1 (Modified Weibull Distribution) 

A random variable X will be said to follow a modified Weibull distribution with exponent c and scale 
parameter x, denoted in the sequel X ~ W(c, x)> if and only if the random variable 



follows a Normal distribution. □ 

These so-called modified-Weibull distributions can be seen to be general forms of the extreme tails of prod- 
uct of random variables (Frisch and Sornette 1997), and using the theorem of change of variable, we can 
assert that the density of such distributions is 



where c and x are the two key parameters. 

These expressions are close to the Weibull distribution, with the addition of a power law prefactor to the 
exponential such that the Gaussian law is retrieved for c = 2. Following Sornette et al. (2000b), Sornette 
et al. (2000a) and Andersen and Sornette (2001), we call (Q the modified Weibull distribution. For c < 1, 
the pdf is a stretched exponential, which belongs to the class of sub-exponential. The exponent c determines 
the shape of the distribution, fatter than an exponential if c < 1. The parameter x controls the scale or 
characteristic width of the distribution. It plays a role analogous to the standard deviation of the Gaussian 
law. 

The interest of these family of distributions for financial purposes have also been recently underlined by 
Brummelhuis and Guegan (2000) and Brummelhuis et al. (2002). Indeed these authors have shown that 




(6) 




(V) 
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given a series of return {r t }t following a GARCH(1,1) process, the large deviations of the returns rt+k 
and of the aggregated returns r t + • • • + r t +k conditional on the return at time t are distributed according 
to a modified- Weibull distribution, where the exponent c is related to the number of step forward k by the 
formula c = 2/k . 

A more general parameterization taking into account a possible asymmetry between negative and positive 
values (thus leading to possible non-zero mean) is 

if x > (8) 
if x < . (9) 



p(x) 
p(x) 



1 c + 

V5F Y TT- 

1 c_ 



d 2 L e Vx+, 



In what follows, we will assume that the marginal probability distributions of returns follow modified 
Weibull distributions. Figurefjshows the (negative) "Gaussianized" returns Y defined in © of the Standard 
and Poor's 500 index versus the raw returns X over the time interval from January 03, 1995 to December 
29, 2000. With such a representation, the modified-Weibull distributions are qualified by a power law of 
exponent c/2, by definition^ The double logarithmic scales of figureQclearly shows a straight line over an 
extended range of data, qualifying a power law relationship. An accurate determination of the parameters 
(X,c) can be performed by maximum likelihood estimation (Sornette 2000, pp 160-162). However, note 
that, in the tail, the six most extreme points significantly deviate from the modified-Weibull description. 
Such an anomalous behavior of the most extreme returns can be probably be associated with the notion 
of "outliers" introduced by Johansen and Sornette (1998, 2002) and associated with behavioral and crowd 
phenomena during turbulent market phases. 

The modified Weibull distributions defined here are of interest for financial purposes and specifically for 
portfolio and risk management, since they offer a flexible parametric representation of asset returns dis- 
tribution either in a conditional or an unconditional framework, depending on the standpoint prefered by 
manager. The rest of the paper uses this family of distributions. 



1.2 Tail equivalence for distribution functions 

An interesting feature of the modified Weibull distributions, as we will see in the next section, is to enjoy the 
property of asymptotic stability. Asymptotic stability means that, in the regime of large deviations, a sum 
of independent and identically distributed modified Weibull variables follows the same modified Weibull 
distribution, up to a rescaling. 

Definition 2 (Tail equivalence) 

Let X and Y be two random variables with distribution function F and G respectively. 

X and Y are said to be equivalent in the upper tail if and only if there exists A + G (0, oo) such that 

lim 1 —^f\ = X + . (10) 

Similarly, X and Y are said equivalent in the lower tail if and only if there exists A_ G (0, oo) such that 

lim ^4 = A_. (11) 

G{x) 



-oo 



□ 
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Applying l'Hospital's rule, this gives immediately the following corollary: 
Corollary 1 

LetX and Y be two random variables with densities functions f and g respectively. X and Y are equivalent 
in the upper (lower) tail if and only if 

Hm 44 = A ±' A±€(0,oo). (12) 

x^±oo g(x) 

□ 

1.3 The Gaussian copula 

We recall only the basic properties about copulas and refer the interested reader to (Nelsen 1998), for in- 
stance, for more information. Let us first give the definition of a copula of n random variables. 

Definition 3 (Copula) 

A function C : [0, l] n — ► [0, 1] is a n-copula if it enjoys the following properties : 

• Vug [0,1], C(l,--- ,1,«,1--- ,1) = u, 

• Viij G [0, 1], C(ui, ■ ■ • , u n ) = if at least one of the Ui equals zero , 

• C is grounded and n-increasing, i.e., the C-volume of every boxes whose vertices lie in [0, 1]™ is 
positive. □ 

The fact that such copulas can be very useful for representing multivariate distributions with arbitrary 
marginals is seen from the following result. 

Theorem 1 (Sklar's Theorem) 

Given an n-dimensional distribution function F with continuous marginal distributions F\, ■ ■ ■ ,F n , there 
exists a unique n-copula C : [0, l] n — > [0, 1] such that : 

F(x lr -- ,x n )=C(F 1 (x 1 ),--- ,F n (x n )) . (13) 

□ 



This theorem provides both a parameterization of multivariate distributions and a construction scheme for 
copulas. Indeed, given a multivariate distribution F with margins F\, ■ ■ ■ , F n , the function 

C(ui,--- ,u n ) = F(F l - 1 (u 1 ),--- ,F-\u n )) (14) 

is automatically a n-copula. Applying this theorem to the multivariate Gaussian distribution, we can derive 
the so-called Gaussian copula. 

Definition 4 (Gaussian copula) 

Let <I> denote the standard Normal distribution and ^v,n the n-dimensional Gaussian distribution with cor- 
relation matrix V. Then, the Gaussian n-copula with correlation matrix V is 

CV(m, ■••,«„) = $v> • • • , $ -1 (wn)) , (15) 

whose density 

dC v (ui,--- ,u n ) 

c v (ui,--- ,u n ) = — (16) 

OUi ■ ■ ■ ou n 
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reads 

c v (u ir -- ,u n ) = -^=L=exp f-^2/( u )(V -1 -Id)y (u) J (17) 

with yk(u) = <I> _1 (ufc). Note that theorem 1 and equation fi4l) ensure that Cy(ui, • • • , u n ) in equation f/5l ) 
is a copula. □ 

It can be shown that the Gaussian copula naturally arises when one tries to determine the dependence be- 
tween random variables using the principle of entropy maximization (Rao 1973, Sornette et al. 2000b, for 
instance). Its pertinence and limitations for modeling the dependence between assets returns has been tested 
by Malevergne and Sornette (2001), who show that in most cases, this description of the dependence can 
be considered satisfying, specially for stocks, provided that one does not consider too extreme realizations 
(Malevergne and Sornette 2002a, Malevergne and Sornette 2002b, Mashal and Zeevi 2002). 



2 Portfolio wealth distribution for several dependence structures 
2.1 Portfolio wealth distribution for independent assets 

Let us first consider the case of a portfolio made of independent assets. This limiting (and unrealistic) case 
is a mathematical idealization which provides a first natural benchmark of the class of portolio return dis- 
tributions to be expected. Moreover, it is generally the only case for which the calculations are analytically 
tractable. For such independepent assets distributed with the modified Weibull distributions, the following 
results prove the asymptotic stability of this set of distributions: 

Theorem 2 (Tail equivalence for led modified Weibull random variables) 

Let Xi, X2, ■ ■ ■ , Xn be N independent and identically W(c, x) -distributed random variables. Then, the 

variable 

S N = X l +X 2 + --- + X N (18) 
is equivalent in the lower and upper tail to Z ~ W(c, x), with 

X = N^ X , c>l, (19) 
X = X, c<l. (20) 

□ 

This theorem is a direct consequence of the theorem stated below and is based on the result given by Frisch 
and Sornette (1997) for c > 1 and on general properties of sub-exponential distributions when c < 1. 

Theorem 3 (Tail equivalence for weighted sums of independent variables) 

Let X\,X2,--- ,Xn be N independent and identically W(c,x) -distributed random variables. Let w\, 

W2 , ■ ■ ■ , wn be N non-random real coefficients. Then, the variable 

Sn = w\X\ + W2X2 H ^wnX n (21) 

is equivalent in the upper and the lower tail to Z ~ W(c, x) with 

X = (XJKI^M -X, Ol, (22) 
X = max{|«;i|, \w 2 \, ■ ■ ■ , \w N \}, c < 1. (23) 

□ 
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The proof of this theorem is given in appendix 1X1 



Corollary 2 

Let Xx, X 2 , ■ ■ ■ , Xjy be N independent random variables such that X,- L ~ W(c, Xi)- Let wx, w 2 , ■ ■ ■ , 
be N non-random real coefficients. Then, the variable 

Sn = wxXx + w 2 X 2 H 1- wnX n (24) 

is equivalent in the upper and the lower tail to Z ~ W(c, x) with 

X = (Yj\ w iXi\^\ , c>l, (25) 
X = max{\w 1 xx\,\w 2 X2\,--- ,\wnXn\}, c < 1. (26) 

□ 

The proof of the corollary is a straightforward application of theorem [5] Indeed, let Yi,Y 2 , 
independent and identically W(c, 1) -distributed random variables. Then, 

(Xx,X 2 ,--- ,X N ) = {xiY\,X2Y 2 ,--- ,XnYn), 

which yields 

Sn = mxi • Yx + w 2 \2 • Y 2 H h watXw ■ Y N . 

Thus, applying theorem|51to the i.i.d variables Yj's with weights WiXi leads to corollary El 

2.2 Portfolio wealth distribution for comonotonic assets 

The case of comonotonic assets is of interest as the limiting case of the strongest possible dependence 
between random variables. By definition, 

Definition 5 (Comonotonicity) 

the variables Xx, X 2 , ■ ■ ■ , Xn are comonotonic if and only if there exits a random variable U and non- 
decreasing functions fx, f2, • • ■ , In such that 

(Xx,X 2 ,--- ,X N )±(fx(U)J 2 (U) r -- J N (U)). (29) 

□ 



• • ■ , Y N be N 
(27) 
(28) 



In terms of copulas, the comonotonicity can be expressed by the following form of the copula 

C(ux,u 2 , ■ ■ ■ ,u N ) = min(ui,u 2 , • • • ,u N ) . (30) 

This expression is known as the Frechet-Hoeffding upper bound for copulas (Nelsen 1998, for instance). It 
would be appealing to think that estimating the Value-at-Risk under the comonotonicity assumption could 
provide an upper bound for the Value-at-Risk. However, it turns out to be wrong, due -as we shall see in the 
sequel- to the lack of coherence (in the sense of Artzner et al. (1999)) of the Value-at-Risk, in the general 
case. Notwithstanding, an upper and lower bound can always be derived for the Value-at-Risk (Embrechts 
et al. 2002b). But in the present situation, where we are only interested in the class of modified Weibull 
distributions with a Gaussian copula, the VaR derived under the comonoticity assumption will actually 
represent the upper bound (at least for the VaR calculated at sufficiently hight confidence levels). 
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Theorem 4 (Tail equivalence for a sum of comonotonic random variables) 

Let Xi,X 2 , ■ ■ ■ , Xjy be N comonotonic random variables such that ~ W(c, Xi)- Let wi,W2, ■ ■ ■ , 

be N non-random real coefficients. Then, the variable 

Sn = wiXi + W2X2 H \-wnX n (31) 

is equivalent in the upper and the lower tail to Z ~ W(c, x) with 

X = J2 WiXi ■ (32) 

i 

□ 

The proof is obvious since, under the assumption of comonotonicity, the portfolio wealth S is given by 

N 



d 

i i=l 

and for modified Weibull distributions, we have 



S = ^2w i -X i ^^2w i -f i (U), (33) 



fi(-) = sgn(-) Xl (^=) ' 1 , (34) 

in the symmetric case while U is a Gaussian random variable. If, in addition, we assume that all assets have 
the same exponent Cj = c, it is clear that S ~ W(c, x) with 



X = ^ w - Xi . (35) 



It is important to note that this relation is exact and not asymptotic as in the case of independent variables. 

When the exponents Cj's are different from an asset to another, a similar result holds, since we can still write 
the inverse cumulative function of S as 

N 

f sHp) = /Z ™i F x-(p), P e (0, 1), (36) 
i=i 

which is the property of additive comonotonicity of the Value-at-Risk 1 . Let us then sort the Xj 's such that 
c\ = C2 = • • • = c p < c p+ i < ■ ■ ■ < c^. We immediately obtain that S is equivalent in the tail to 
Z ~ W(ci, x). where 

p 

X = J2wiXi- (37) 

i=i 

In such a case, only the assets with the fatest tails contributes to the behavior of the sum in the large deviation 
regime. 



'This relation shows that, in general, the VaR calculated for comonotonic assets does not provide an upper bound of the VaR, 
whatever the dependence structure the portfolio may be. Indeed, in such a case, we have VaR(Xi + X2) = VaR(X 1 ) + VaR^A^) 
while, by lack of coherence, we may have VaR(Ai + X2) > VaR(Ai) + VaR(A2) for some dependence structure between X\ 
and X2 ■ 
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2.3 Portfolio wealth under the Gaussian copula hypothesis 



2.3.1 Derivation of the multivariate distribution with a Gaussian copula and modified Weibull mar- 
gins 

An advantage of the class of modified Weibull distributions © is that the transformation into a Gaussian, 
and thus the calculation of the vector y introduced in definition [0 is particularly simple. It takes the form 



y k = sgn(xfc) y/2 ( — 

Xk 



(38) 



where is normally distributed . These variables Yi then allow us to obtain the covariance matrix V of the 
Gaussian copula : 



Vij = 2 • E 



1 Xi 



T ■ I \ 2 / 7* • \ 2 

tlLt'l, | \ / | •X' J | 



X,j 



(39) 



which always exists and can be efficiently estimated. The multivariate density P(x) is thus given by: 



A' 



P(xi,--- ,x N ) = c v (xi,x 2 ,--- ,%n) IJpi(^i) 



i=l 



1 w 

2V/ 2 v^n 

v i=i 



c% x% 



.|c/2-l 



c/2 



■ exp 



c/2 



c/2 



(40) 
(41) 



Obviously, similar transforms hold, mutatis mutandis, for the asymmetric case d8l9t . 



2.3.2 Asymptotic distribution of a sum of modified Weibull variables with the same exponent c > 1 

We now consider a portfolio made of dependent assets with pdf given by equation ATI or its asymmetric 
generalization. For such distributions of asset returns, we obtain the following result 

Theorem 5 (Tail equivalence for a sum of dependent random variables) 
Let Xi, X2, ■ ■ ■ , Xn be N random variables with a dependence structure described by the Gaussian copula 
with correlation matrix V and such that each X{ ~ W(c,Xi)- Let wi, W2, ■ ■ ■ , wn be N (positive) non- 
random real coefficients. Then, the variable 



Sn = W1X1 + W2X2 H 1- wnX n (42) 

is equivalent in the upper and the lower tail to Z ~ W(c, x) with 

c— 1 

where the cj{ 's are the unique (positive) solution of 

J2 V ik 1(J i C/2 = w kXk °\7 C ' 2 , Vfc . (44) 

i 

□ 
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X 


A_ 


Independent Assets 


c-l 

(EiliKxil^ 1 ) c ,o i 

max{|wixi|, • • • , \w n xn\}, c < 1 


AT-l 

r c i 2 

Card {\wiXi\ = maxj{ \wjXj\}} 


Comonotonic Assets 




1 


Gaussian copula 


c-l 

(Ei^iXiO-i) c , c > 1 


see appendix IE1 



Table 1: Summary of the various scale factors obtained for different distribution of asset returns. 

The proof of this theorem follows the same lines as the proof of fheorem[3] We thus only provide a heuristic 
derivation of this result in appendix |B] Equation (l44l > is equivalent to 

v ik WkXk Ok 1 ^' 2 = of' 2 , Vi . (45) 

k 

which seems more attr active since it does not require the inversion of the correlation matrix. In the special 

case where V is the identity matrix, the variables Xj's are independent so that equation d43t must yield the 

i 

same result as equation d22l . This results from the expression of = (wkXk) c_1 valid in the independent 

case. Moreover, in the limit where all entries of V equal one, we retrieve the case of comonotonic assets. 

Obviously, V -1 does not exist for comonotonic assets and the derivation given in appendix|B]does not hold, 

i 

but equation (I45t remains well-defined and still has a unique solution = w^Xk) c_1 which yields the 

scale factor given in theorem|4] 
2.4 Summary 

In the previous sections, we have shown that the wealth distribution F${x) of a portfolio made of assets with 
modified Weibull distributions with the same exponent c remains equivalent in the tail to a modified Weibull 
distribution W(c,x). Specifically, 

F s {x) ~ A_ F z {x) , (46) 

when x — > — oo, and where Z ~ W(c, x)- Expression d46l defines the proportionality factor or weight A_ 
of the negative tail of the portfolio wealth distribution Fg(x). Table Q summarizes the value of the scale 
parameter x f° r the different types of dependence we have studied. In addition, we give the value of the 
coefficient A_, which may also depend on the weights of the assets in the portfolio in the case of dependent 
assets. 
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3 Value-at-Risk 



3.1 Calculation of the VaR 

We consider a portfolio made of N assets with all the same exponent c and scale parameters Xi> i £ 
{1, 2, • • • , N}. The weight of the i th asset in the portfolio is denoted by Wi. By definition, the Value- 
at-Risk at the loss probability a, denoted by VaR a , is given , for a continuous distribution of profit and loss, 
by 

Pr{W(r) - W{0) < -VaR a } = a, (47) 

which can be rewritten as 

VaR a I 



Pr i s< ^} = "- <48) 

In this expression, we have assumed that all the wealth is invested in risky assets and that the risk-free 
interest rate equals zero, but it is easy to reintroduce it, if necessary. It just leads to discount VaR a by the 
discount factor 1/(1 + /j,q), where jjlq denotes the risk-free interest rate. 

Now, using the fact that F${x) ~ A_ Fz(x), when x — > — oo, and where Z ~ W(c, x), we have 



as VaR a goes to infinity, which allows us to obtain a closed expression for the asymptotic Value-at-Risk 
with a loss probability a: 

2/c 



VaR a ~ W(0) V 



1 / OL 

r,1 -r 



(50) 



2Vc 

~ t(afl c W(0) ■ x, (51) 
where the function $(•) denotes the cumulative Normal distribution function and 

«„) , 1,- (l - «.) . (52) 

In the case where a fraction wo of the total wealth is invested in the risk-free asset with interest rate /j,q, the 
previous equation simply becomes 

VaR a ~ i{afl c (1 - w ) ■ W(0) ■ x ~ ^0^(0)^0- (53) 

Due to the convexity of the scale parameter % the VaR is itself convex and therefore sub-additive. Thus, for 
this set of distributions, the VaR becomes coherent when the considered quantiles are sufficiently small. 

The Expected-Shortfall ES a , which gives the average loss beyond the VaR at probability level a, is also 
very easily computable: 

ES a = - ! VaR u (iu (54) 
a Jo 

= ((a)(l-w )-W(0)-x-w W(0)M, (55) 

where ((a) = ^ J Q a £(u) 2 / c du . Thus, the Value-at-Risk, the Expected-Shortfall and in fact any downside 
risk measure involving only the far tail of the distribution of returns are entirely controlled by the scale 
parameter x- We see that our set of multivariate modified Weibull distributions enjoy, in the tail, exactly the 
same properties as the Gaussian distributions, for which, all the risk measures are controlled by the standard 
deviation. 
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3.2 Typical recurrence time of large losses 



Let us translate these formulas in intuitive form. For this, we define a Value-at-Risk VaR* which is such that 
its typical frequency is I/Tq. Tq is by definition the typical recurrence time of a loss larger than VaR*. In 
our present example, we take To equals 1 year for example, i.e., VaR* is the typical annual shock or crash. 
Expression d?9l then allows us to predict the recurrence time T of a loss of amplitude VaR equal to j3 times 
this reference value VaR* : 



Figure |2] shows In jr versus f3. Observe that T increases all the more slowly with j3, the smaller is the 
exponent c. This quantifies our expectation that large losses occur more frequently for the "wilder" sub- 
exponential distributions than for super-exponential ones. 

4 Optimal portfolios 

In this section, we present our results on the problem of the efficient portfolio allocation for asset distributed 
according to modified Weibull distributions with the different dependence structures studied in the previous 
sections. We focus on the case when all asset modified Weibull distributions have the same exponent c, as 
it provides the richest and more varied situation. When this is not the case and the assets have different 
exponents Cj, i = 1, N, the asymptotic tail of the portfolio return distribution is dominated by the asset 
with the heaviest tail. The largest risks of the portfolio are thus controlled by the single most risky asset 
characterized by the smallest exponent c. Such extreme risk cannot be diversified away. In such a case, for 
a risk-averse investor, the best strategy focused on minimizing the extreme risks consists in holding only the 
asset with the thinnest tail, i.e., with the largest exponent c. 

4.1 Portfolios with minimum risk 

Let us consider first the problem of finding the composition of the portfolio with minimum risks, where 
the risks are measured by the Value-at-Risk. We consider that short sales are not allowed, that the risk free 
interest rate equals zero and that all the wealth is invested in stocks. This last condition is indeed the only 
interesting one since allowing to invest in a risk-free asset would automatically give the trivial solution in 
which the minimum risk portfolio is completely invested in the risk-free asset. 

The problem to solve reads: 



In some cases (see table0, the prefactor defined in d52t also depends on the weight Wj's through A_ 
defined in d46t . But, its contribution remains subdominant for the large losses. This allows to restrict the 
minimization to x instead of £(a) 2//c • X- 




(56) 



VaR*, = minVaR a = £(a) 2/c W(0) ■ 

z2 i= i wi = i 



minx 



(57) 
(58) 
(59) 



Wi > Mi. 
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4. 1. 1 Case of independent assets 



"Super-exponential" portfolio (c > 1) 

Consider assets distributed according to modified Weibull distributions with the same exponent c > 1. 
The Value-at-Risk is given by 



c-l 

N \ ~ 



VaR Q = i(af c W(0) ■ £ \w iXi \^ , (60) 



\j=i 

Introducing the Lagrange multiplier A, the first order condition yields 

dx A 



dwi f (a) W(0) 
and the composition of the minimal risk portfolio is 



Vi, (61) 



Wi * = -^— c (62) 



9 2 X 



which satistifies the positivity of the Hessian matrix Hjk — Q w g Wk 
The minimal risk portfolio is such that 



(second order condition). 



V^ = «2™, ,63) 

(E, x7 c ) ' ^ Xj 

where /i, is the return of asset z and ^* is the return of the minimum risk portfolio. 



sub-exponential portfolio (c < 1) 

Consider assets distributed according to modified Weibull distributions with the same exponent c < 1. 
The Value-at-Risk is now given by 

VaR a = £(a) c/2 W(0) • max{\w lX i\, ■ ■ ■ , \w nX n\}- (64) 

Since the weights Wi are positive, the modulus appearing in the argument of the max() function can be 
removed. It is easy to see that the minimum of VaR Q is obtained when all the WiXi's are equal, provided 
that the constraint ^ iUj = 1 can be satisfied. Indeed, let us start with the situation where 

wixi = W2X2 = • • • = wnXn ■ (65) 

Let us decrease the weight w\. Then, w\X\ decreases with respect to the initial maximum situation d63t but, 
in order to satisfy the constraint J2i w i = 1> at l east one 01 tne other weights Wj, j > 2 has to increase, so 
that WjXj increases, leading to a maximum for the set of the WiXi's greater than in the initial situation where 
d65l holds. Therefore, 

w* = — , Vi, (66) 

Xi 

and the constraint ^ Wi = 1 yields 

^4 = ^-T, (6V) 
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and finally 



VaR„ - 'J,"" . // - — ■ (68. 



£(«) c/2 W(0) , = ZiXT 1 * 

E-i ■ Q -1 ' r 1 v-^ — 

j / j% X% / jj Xj 



The composition of the optimal portfolio is continuous in c at the value c = 1. This is the consequence 
of the continuity as a function of c at c = 1 of the scale factor \ f° r a sum °f independent variables. In 
this regime c < 1, the Value-at-Risk increases as c decreases only through its dependence on the pref actor 
£(a) 2//c since the scale factor X remains constant. 

4.1.2 Case of comonotonic assets 

For comonotonic assets, the Value-at-Risk is 

VaR a = £(a) c / 2 W{V) • ]T w iX i (69) 

i 

which leads to a very simple linear optimization problem. Indeed, denoting \i = m i n {xi> Xii"' > Xn}, 
we have 

w iXi >Xi ^2wi = xi, (70) 

i i 

which proves that the composition of the optimal portfolio is w* = 1, w* = i > 2 leading to 

VaR; = e(a) c / 2 iy(0) X i, /i* = Ml- (71) 

This result is not surprising since all assets move together. Thus, the portfolio with minimum Value-at-Risk 
is obtained when only the less risky asset, i.e., with the smallest scale factor Xu is held. In the case where 
there is a degeneracy in the smallest x °f order p (xi = X2 = ••• = Xp = mm {xi;X2> • ■ ■ ,Xn}), the 
optimal choice lead to invest all the wealth in the asset with the larger expected return fij, j G {1, • ■ ■ ,p}. 
However, in an efficient market with rational agents, such an opportunity should not exist since the same 
risk embodied by X i = X2 = ••• = Xp should be remunerated by the same return = = ... = /x p . 



4.1.3 Case of assets with a Gaussian copula 

In this situation, we cannot solve the problem analytically. We can only assert that the miminization problem 
has a unique solution, since the function NaR a {{wi}) is convex. In order to obtain the composition of the 
optimal portfolio, we need to perform the following numerical analysis. 

It is first needed to solve the set of equations J2i V-ij 1 ^ 2 = w jXj a ] ^ or tne equivalent set of equations 
given by <45b . which can be performed by Newton's algorithm. Then one have the minimize the quantity 
Yl w iXi< J i{{ w i})- To this aim, one can use the gradient algorithm, which requires the calculation of the 
derivatives of the crj's with respect to the w^s. These quantities are easily obtained by solving the linear set 
of equations 

^ d^ + {-2- l ) w ^v J d^ = x ^^ (72) 

i J 
Then, the analytical solution for independent assets or comonotonic assets can be used to initialize the 
minimization algorithm with respect to the weights of the assets in the portfolio. 
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4.2 VaR-efficient portfolios 



We are now interested in portfolios with minimum Value-at-risk, but with a given expected return \i = 
uij/Xj. We will first consider the case where, as previously, all the wealth is invested in risky assets and 
we then will discuss the consequences of the introduction of a risk-free asset in the portfolio. 



4.2.1 Portfolios without risky asset 

When the investors have to select risky assets only, they have to solve the following minimization problem: 

VaR* = minVaR a = £(a) W(0) ■ mm X (73) 
Eill w i^i = M (74) 

££i m = 1 (75) 
Wi > V*. (76) 

In contrast with the research of the minimum risk portfolios where analytical results have been derived, we 
need here to use numerical methods in every situation. In the case of super-exponential portfolios, with or 
without dependence between assets, the gradient method provides a fast and easy to implement algorithm, 
while for sub-exponential portfolios or portfolios made of comonotonic assets, one has to use the simplex 
method since the minimization problem is then linear. 

Thus, although not as convenient to handle as analytical results, these optimization problems remain easy to 
manage and fast to compute even for large portfolios. 



4.2.2 Portfolios with risky asset 

When a risk-free asset is introduced in the portfolio, the expression of the Value-at-Risk is given by equation 
(1531 . the minimization problem becomes 

VaR* = min£(a) 2 / c (1 - w ) ■ W(0) ■ x ~ w W{0)fi (77) 

Ej=i mvi = v (78) 

£,=i^ = l (79) 
Wi > Vi (80) 

When the risk-free interest rate /j,q is non zero, we have to use the same numerical methods as above to 
solve the problem. However, if we assume that = 0, the problem becomes amenable analytically. Its 
Lagrangian reads 

£ = £(a) 2/c (l-™ )-W(0)-x-Ai lj2w ifH -fj\ - A 2 (^^-lV (81) 

= £(«) 2/c l^wA -W^-x-Xi l^Wi^-A, (82) 
\j¥*> ) \¥o J 

which allows us to show that the weights of the optimal portfolio are 

* /I ^ A u n* ^(«) 2/C (1 ~ W ) ■ W(0) 

Wi = (1 - w ) ■ —jf and VaR Q = (83) 

E,=i wj 2 
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where the wis are solution of the set of equations 




Expression (l%3t shows that the efficient frontier is simply a straight line and that any efficient portfolio is 
the sum of two portfolios: a "riskless portfolio" in which a fraction wo of the initial wealth is invested and 
a portfolio with the remaining (1 — wq) of the initial wealth invested in risky assets. This provides another 
example of the two funds separation theorem. A CAPM then holds, since equation dS^T l together with the 
market equilibrium assumption yields the proportionality between any stock return and the market return. 
However, these three properties are rigorously established only for a zero risk-free interest rate and may not 
remain necessarily true as soon as the risk-free interest rate becomes non zero. 

Finally, for practical purpose, the set of weights w*'s obtained under the assumption of zero risk-free interest 
rate jXQ, can be used to initialize the optimization algorithms when jiQ does not vanish. 



5 Conclusion 



The aim of this work has been to show that the key properties of Gaussian asset distributions of stability 
under convolution, of the equivalence between all down-side riks measures, of coherence and of simple 
use also hold for a general family of distributions embodying both sub-exponential and super-exponential 
behaviors, when restricted to their tail. We then used these results to compute the Value-at-Risk (VaR) and 
to obtain efficient porfolios in the risk-return sense, where the risk is characterized by the Value-at-Risk. 
Specifically, we have studied a family of modified Weibull distributions to parameterize the marginal distri- 
butions of asset returns, extended to their multivariate distribution with Gaussian copulas. The relevance to 
finance of the family of modified Weibull distributions has been proved in both a context of conditional and 
unconditional portfolio management. We have derived exact formulas for the tails of the distribution P(S) 
of returns S of a portfolio of arbitrary composition of these assets. We find that the tail of P(S) is also 
asymptotically a modified Weibull distribution with a characteristic scale x function of the asset weights 
with different functional forms depending on the super- or sub-exponential behavior of the marginals and 
on the strength of the dependence between the assets. The derivation of the portfolio distribution has shown 
the asymptotic stability of this family of distribution with the important economic consequence that any 
down-side risk measure based upon the tail of the asset returns distribution are equivalent, in so far as they 
all depends on the scale factor x an d keep the same functional form whatever the number of assets in the 
portfolio may be. Our analytical study of the properties of the VaR has shown the VaR to be coherent. This 
justifies the use of the VaR as a coherent risk measure for the class of modified Weibull distributions and 
ensures that portfolio optimization problems are always well-conditioned even when not fully analytically 
solvable. The Value-at-Risk and the Expected-Shortfall have also been shown to be (asymptotically) equiv- 
alent in this framework. In fine, using the large class of modified Weibull distributions, we have provided 
a simple and fast method for calculating large down-side risks, exemplified by the Value-at-Risk, for assets 
with distributions of returns which fit quite reasonably the empirical distributions. 
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A Proof of theorem |3] : Tail equivalence for weighted sums of modified 
Weibull variables 

A.l Super-exponential case: c > 1 

Let X\,X2, ■ ■ ■ ,X]\f be N i.i.d random variables with density p(-). Let us denote by /(•) and <?(•) two 
positive functions such that p(-) = g(-) ■ e~^'\ Let w\,W2, ■ ■ ■ ,wn be ./V real non-random coefficients, 

and S = J2i=i w i x i- 



(85) 



Let X = {x G M. , Ysi=i w i x i = The density of the variable S is given by 

P S (S) = [ rfxe^w^- 1115 ^'!, 
Jx 



We will assume the following conditions on the function / 

1. /(■) is three times continuously differentiable and four times differentiable, 

2. f( 2 '{x) > 0, for |x| large enough, 



3. limbic 



/( 3 )(a) 



o, 



4. /^ 3 ) is asymptotically monotonous, 

5. there is a constant (3 > 1 such that ^J^r remains bounded as x goes to infinity, 

6. g(-) is ultimately a monotonous function, regularly varying at infinity with indice v. 

Let us start with the demonstration of several propositions. 



Proposition 1 

under hypothesis 3, we have 

□ 



lim \x\ • f"(x) = 0. 

x— >±oo 



(86) 



Proof 



Hypothesis 3 can be rewritten as lim 



d 1 



x^ztoo dx /( 2 ) {x) 
Ve > 0,3A £ /x > A £ 



0, so that 
d 1 



dx f( 2 )(x) 



< e. 



(87) 



Now, since f" is differentiable, 1 / /" is also differentiable, and by the mean value theorem, we have 

1 1 



f"(x) f"(y) 





d 1 


= \x-y\- 


dU"(0 



(88) 



for some ^ G (x, y). 
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Choosing x > y > A e , and applying equation d87l together with d88l yields 



f"(x) f»(y) 
Now, dividing by x and letting x go to infinity gives 



which concludes the proof. □ 

Proposition 2 

Under assumption 3, we have 



□ 



lim 

x— >oo 



1 



X ■ f"(x) 



< e ■ \x — y\. 



<e, 



lim / (x) = +oo. 

x— >±oo 



(89) 



(90) 



(91) 



Proof 

According to assumption 3 and proposition lim x • f"{x) = do, which means 

a;— »±oo 



Va > 0, 3^ Q /x > A e => x • /"(a?) > a. 



This thus gives 



Vx > a a , x • f"(x) > a 



> - 

X 

f' x f' x dt 

/ f"(t)dt>a. j 

/'(a?) > a • In x — a • In A a + /'(Ax). 



(92) 

(93) 
(94) 
(95) 



The right-hand-side of this last equation goes to infinity as x goes to infinity, which concludes the proof. □ 



Proposition 3 

Under assumptions 3 and 6, the function g(-) satisfies 



C 



uniformly in h, for any positive constant C. □ 
Proof For g non-decreasing, we have 

c a(x(i 



v 9{x + h) 
x^±oo q[x) 



(96) 



c 



x-f"(x) 



< g{x + h) < 9 _X C \} + W%) 



(97) 



f"(x) g(x) g{x) g(x) 

If g is non-increasing, the same inequalities hold with the left and right terms exchanged. Therefore, the 
final conclusion is easily shown to be independent of the monotocity property of g. From assumption 3 and 
proposition H we have 

Va > 0, 3A a /x > A e =^ x ■ f"(x) > a. (98) 
Thus, for all x larger than A a and all \h\ < C/f"(x) 



g (*(!-§)) < 9± + V < 9(x(l + %)) 



g(x) 



g{x) 



g(x) 



(99) 
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Now, letting x go to infinity, 



!_£)"< Km ^<( 1+ ^y, (loo, 



for all a as large as we want, which concludes the proof. □ 
Proposition 4 

Under assumptions 1, 3 and 4 we have, for any positive constant C: 



wm^ C v s m^,x+h]\f {3) (0\ 

y\h\ < -Tiff^, ton -—2 = 0. (101) 



□ 



Proof 

Let us first remark that 



su P^e[x,x+h] 


/ (3) (0| 




/(3)(x) 





< 



|/ (3 >N 



/(3) (x 


(} + x/"(r)) 


) 




/(3)(x) 





/ (3) (^(1+g)) 

|/(3)(x)| 



(103) 

(104) 
(105) 



|/ (3) (e)| 

f»{x)* |/(3)(x)| ' /"(^) 2 ' ^ ^ 

The rightmost factor in the right-hand-side of the equation above goes to zero as x goes to infinity by as- 
sumption 3. Therefore, we just have to show that the leftmost factor in the right-hand-side remains bounded 
as x goes to infinity to prove Proposition 4. 

Applying assumption 4 according to which /( 3 ) is asymptotically monotonous, we have 



for every x larger than some positive constant A a by assumption 3 and proposition ^ Now, for a large 
enough, 1 + ^ is less than [3 (assumption 5) which shows that ~~^^n^)j — remains bounded for large 
x, which conclude the proof. □ 

We can now show that under the assumptions stated above, the leading order expansion of Ps{S) for large 
S and finite N > 1 is obtained by a generalization of Laplace's method which here amounts to remark that 
the set of x* 's that maximize the integrand in d85t are solution of 

f&xt) = a(S) Wi , (106) 

where cr(S) is nothing but a Lagrange multiplier introduced to minimize the expression Y2i=i fi( x i) under 
the constraint J2iLi w -i x i = S. This constraint shows that at least one Xi, for instance x\, goes to infinity 
as S — ► oo. Since /'(a?i) is an increasing function by assumption 2 which goes to infinity as x\ — * +oo 
(proposition expression (1106b shows that cr(S) goes to infinity with S, as long as the weight of the asset 
1 is not zero. Putting the divergence of a(S) with S in expression (1106b for i = 2, N ensures that each 
x* increases when S increases and goes to infinity when S goes to infinity. 
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Expanding fi(xi) around x* yields 



rx*+hi ft 

f{xi) = f{xt) + f(x*) -hi+ I dt du /"(«) 



(107) 



where the set of hi = X{ — x* obey the condition 



N 



Wihi = . 



(108) 



i=l 



Summing dl06t in the presence of relation d 108l> . we obtain 

N N rxf+hi 



£ f(xi) = £ f(x*) +J2 dt du f"(u) . 

i=\ i=i i=i Jx t Jx * 



(109) 



Thus exp(— ^2f{xi)) can be rewritten as follows : 



exp 



N 



8=1 



exp 



N N f X *+h- ft 



i=l 



i=l 



(110) 



Let us now define the compact set Ac = {h G W N , X)i=i f"( x i ) 2 ' — C 2 } f° r an y gi ven positive 
constant C and the set H = {hG M. N , ^ i=1 = 0}. We can thus write 

P s (5) = [ dhe-^ [fix ^- ln9 ^ )] , 
Jh 

= [ dhe-^M^- 1119 ^ + [ dhe-^^ 1 ^- 1119 ^ , 
JAcnn JAcnH 

We are now going to analyze in turn these two terms in the right-hand-side of ( 11121 . 

First term of the right-hand-side of d 1 1 2b . 

Let us start with the first term. We are going to show that 



(111) 
(112) 



lim 

S— »oo 



■ Ell lT hl dtj*. du /"(«)-ln fl (xi) 



N-l 

(2tt)^- 11,3(0 

W 2 JJ^ 1 /"(**) 
N I J'f 1 3 3 



1, for some positive C. (113) 



In order to prove this assertion, we will first consider the leftmost factor of the right-hand-side of ( II 12 

-Y.fi +hl */*. du /"(«) 



•E/(*0 



(114) 



-i— r . , . lx-t»(„*\ h 2 -T. T* l+hi dt f\ du[f"(u)-f"(x*)] 

= Yl9(x* + hi) e~2 E/(»i)Ai e J ^ v . (115) 

Since for all (Gl, e - ^ < < e^l, we have 

E/,7'''''" 1 1: du [/"(«)-/"(xJ 



, -E O + " 1 d * /-* du [/"(«)-/"(*<)] ^ 
< e i 1 < e 



ZO +h *dtf^ du [/'»-/"«)] 

i i 

(116) 
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and 



-ZJ*i +hi dtf**du\f"(u)-f"(x*)\ -ZO^dtJ* t dulf"(u)-f"(x*)] ZO^dt&du\f"(u)-f"(x*)\ 

e » 1 <. e * 1 <. e * 1 , 

(117) 

since whatever the sign of hi, the quantity G 1 dt f f <iu — /"(^*)| remains always positive. 

But, \u — a?|| < \hi\ < jrfj^j, which leads, by the mean value theorem and assumption 1, to 



r x i + h i rt 



l/» -/"(**)! < S up i/( 3 )(e)i 



U — Xi 



C 



which yields 



< SUp l/W (0| jji—,, 

£<=(x*,x*+hi) J \ x i) 

< sup|/( 3 )(e)|-^— , where ft 



0<£/ dt / d W |/»- /''(<)!< sup |/( 3 )(0I 7 ^-—^. (121) 



c c 

1 Xj T" 



/"(*?)' 4 /"«) 



(118) 
(119) 
(120) 



Thus 



I E 4 sup l/W (0| ^fc? -E;.*/ +hl */*.du[/»-/"(x?)] ^ IE, sup 1/(3) (0| jr^-h} 



< e' 



(122) 



where sup^gg. |/^(0I> have been denoted by sup|/( 3 )(0l in the previous expression, in order not to 
cumber the notations. 



By proposition @1 we know that for all h £ Ac and all > 

sup|/( 3 )(£) 



< £j, for x* large enough, 



(123) 



so that 



, C-e' ^ h 2 -E f X l +hi dt I - ** du [/"(«)-/" C-e' ^ h 2 

Ve' > and Vh G „4 C , e - — E * h » < e 3x * » 1 < h * , 

for |x| large enough. 

Moreover, from proposition |3] we have for all £j > and x* large enough: 



(124) 



<0- + <*)", 



so, for all e" > 



VheAc, (l-e>Y u <T\ 9iX ) + * M) <(HeT. 



(125) 



(126) 



Then for all e > and |x| large enough, this yields : 



(1 _ e \Nu e -\ Ei(/"(*?)+o- e )-ft? < TT g( a «) e - E /(**) < (1 + e )JV«/ e -| E 4 (/" W)-c-e)-h? (12 7) 
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for all h G Ac- Thus, integrating over all the h £ Ac H 7i and by continuity of the mapping 

G(Y) = / dh 5 (h,Y) (128) 

JA c nH 

where g(h, Y) = e~5 £ y i- /l f ( we can conclude that, 

Yl9(x*) f AcnH dhe-^f"(*m ' 

Now, we remark that 



1. (129) 



/h JAcnn JAcnH 

with 



(130) 



/ dhe-iS« = (2?r) 2 ~ - (131) 



and 

JA c nH v ' 

where x* = max{x*} (note that 1/ f"(x) — > oo with x by Proposition 1). Indeed, we clearly have 



(132) 



/ dhe-* < [_dhe-^f"^ h l (133) 

JAcnn J Ac 



(2n) N / 2 



i «? 



du ff - , (134) 

11 /r. JJT7Z*\ ' 



where we have performed the change of variable U{ = f"(x*) ■ hi and denoted by Be the set {h G 
R N , uf < C 2 }. Now, let x^ ax = max{x*} and x^ in = min{x*}. Expression (I134t then gives 

' ^u— (135) 



r d he-iE/"(« < / 

tew " /"(^j^ 2 is 



f"«J N/2 J Be- (27r/"K in ))^' 

" Sn -' m~F r U'^cjJ' (136) 

TV i 

~ 5 ■/ l X max; / ^ . -3 JTT^^ (137) 

which decays exponentially fast for large 5 (or large x* nax ) as long as /" goes to zero at infinity, i.e, for any 



function / which goes to infinity not faster than x . So, finally 



dhe-jE /"(«?)*?= ^ 2 + o( e -7TO) , (138) 

which concludes the proof of equation (^Dj>. 



N-l 
2 
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Second term of the right-hand-side of dll2t . 

We now have to show that 

f dhe-^ f ^ +hi) ~ 9 ^ +hi) (139) 
JAcrM. 

can be neglected. This is obvious since, by assumption 2 and 6, the function f(x) — lng(x) remains convex 
for x large enough, which ensures that f(x) — In g(x) > C%\x\ for some positive constant C\ and x large 
enough. Thus, choosing the constant C in Ac large enough, we have 

[ dhe-^f^- 1 " 19 ^ < f_ dhe- Cl Zti\**i+hi\ „ o (V/ 77 !^) ) . (140) 
JAcnn JAcr\H V / 

Thus, for S large enough, the density Ps{S) is asymptotically equal to 



Ps(S)=Ug(x*)^2l, 



JV-l 
2 



(141) 



7TR) 

In the case of the modified Weibull variables, we have 

/(X)= (^) C ' (M2) 

and 

g(s) = ^ ,„ • M^" 1 , (143) 

which satisfy our assumptions if and only if c > 1. In such a case, we obtain 

i 

x* = i-g- • 5, (144) 

which, after some simple algebraic manipulations, yield 

JV-l 

|S|§ _1 e (145) 



2(c-l)J 2/^Fx c/2 



c-l 



* = c (146) 



P(5) 

with 

as announced in theorem |3] 

A.2 Sub-exponential case: c < 1 

Let Xi,X2, • • • , X/v be TV i.i.d sub-exponential modified Weibull random variables W(c, %), with distribu- 
tion function F. Let us denote by Gs the distribution function of the variable 

Sn = wiXi + w 2 X 2 H h wnX n , (147) 

where wi , , • • • , it>7v are real non-random coefficients. 

Let w* = max{\wi\i \ w 2\, ••• , \wn\}. Then, theorem 5.5 (b) of Goldie and Kliippelberg (1998) states that 

lim ° s{ l/ W - = Card {% e {1, 2, , TV} : lioJ = w*}. (148) 
x^oo .r (a:) 

By definitional this allows us to conclude that Sn is equivalent in the upper tail to Z ~ W(c, itf*x). 
A similar calculation yields an analogous result for the lower tail. 



24 



B Asymptotic distribution of the sum of Weibull variables with a Gaussian 
copula. 



We assume that the marginal distributions are given by the modified Weibull distributions: 



P% \X% 



1 C 

V5F^ 



(149) 



and that the Xi' s are all equal to one, in order not to cumber the notation. As in the proof of corollary |2j it 
will be sufficient to replace u>i by WiXi to reintroduce the scale factors. 

Under the Gaussian copula assumption, we obtain the following form for the multivariate distribution : 



P(xi, ■ ■ ■ ,x N ) 



jv 



N 



nc/2-1 
x( exp 



2V/ 2 v / drtF j i = i 



En 



-1 c/2 c/2 



i -J 



(150) 



Let 



c/2 _ .c/2 



(151) 



'•J 



We have to minimize / under the constraint ^ = S. As for the independent case, we introduce a 
Lagrange multiplier A which leads to 



(152) 



The left-hand-side of this equation is a homogeneous function of degree c — 1 in the x*'s, thus necessarily 



(153) 



where the ov's are solution of 



V jk °3 ' a k 1 =W k 



(154) 



The set of equations (I154t has a unique solution due to the convexity of the minimization problem. This 
set of equations can be easily solved by a numerical method like Newton's algorithm. It is convenient to 
simplify the problem and avoid the inversion of the matrix V, by rewritting ( 11541 as 



Vjk w k a k 

k 



l-c/2 = a c/2 



(155) 



Using the constraint ^ WiX* = S, we obtain 



S 



so that 



■ s. 



(156) 



(157) 
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Thus 



f( Xl ,---, XN ) = f(x* ir -- , x * N ) + Y^^-h i + \Y.^rLr h i h o + 



■ij 



s c 



where, as in the previous section, /ij = Xj — x* and the derivatives of / are expressed at x* , 

It is easy to check that the n* A -order derivative of / with respect to the Xj's evaluated at {x*} 
to S c ~~ n . In the sequel, we will use the following notation : 



d n f 



= M (n) . S c - n 



We can write : 

/(xi,--- ,x N ) = 



up to the fourth order. This leads to 

g c T S c_2 (2) / \ 

P(5) oc e C »^,) c - 1 dh!--- dh N e — M ^ hih > 6 (Y^ w i h i) 



S 



c-3 



1 + £ 7-E M S^^ + 



ijk 



Using the relation 5 wA) = / f§ e ifc ^' we obtain : 



S 



c-3 



i + VE M S^* + 



or in vectorial notation 



P(S) oc e C 



ijfc 



dh e -^-- htM(2)h - iA:wth x 



-c-3 



ijk 



Let us perform the following standard change of variables 



h = h' 



ik 
S^2 



w . 



(M^ 2 ) 1 exists since / is assumed convex and thus M^ 2 ) positive) : 



— h*M( 2 )h + ikv^h = — h"M( J )h' + _*_ w t M (2)-i 

2 2 2S' C - 2 
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This yields 



P(S) oc e E ' 



/ 



x dhe 2 



^o 1 ( h+^MirVj'MP) f h+^M< 2 ) 



' s c -' 2 



5" 



c-3 



2vr 



(167) 



Denoting by (-)h the average with respect to the Gaussian distribution of h and by (-) k the average with 
respect to the Gaussian distribution of k, we have : 



P(S) - \l ^ M ^r (2nS^)^ e'^k^ x 



?c-3 



wtrvK 2 ) w 



(4) 



ijk 
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ijkl 



(168) 



We now invoke Wick's theorem 2 , which states that each term ((hi ■ ■ ■ h p )h)k can be expressed as a product 
of pairwise correlation coefficients. Evaluating the average with respect to the symmetric distribution of k, 
it is obvious that odd-order terms will vanish and that the count of powers of S involved in each even-order 
term shows that all are sub-dominant. So, up to the leading order : 



P(S) - \l d T^-T (2^-)^ e -(E^F- . 



w 



tM( 2 )" 



w 



(169) 



The matrix can be calculated, which yields 



M, 



(2) _ 



1 



kl 



and shows that 



' det M( 2 ) 
w t M (2)- i w 



(at-i)(|-i) / detM- 



w t M" 1 w 



(170) 
(171) 

(172) 



The inverse matrix M 1 satisfies Yli M k i ■ (M = 5 k j which can be rewritten: 



(173) 



or equivalently 



2 

(| - l) w k (M- l ) k] + C -Y, Vu X ■ (M-%4 af" 1 = S kj ■ a k (174) 

i 



See for instance (Brezin et al. 1976) for a general introduction, (Sornette 1998) for an early application to the portfolio problem 
and (Sornette et al. 2000b) for a systematic utilization with the help of diagrams. 
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which gives 



c > ~ 2 



1) e ™ k {M-%i w, + \ E ^t 1 • ^ 3 ^ = E ^ • w • (175) 



j,k j,k,l j,k 

Summing the rightmost factor of the left-hand-side over k, and accounting for equation d!54t leads to 

2 

c (J " X ) E Wk (^~% Wj + — E «i {M~ 1 )ijWj = E <W ■ (176) 



so that 

w t M _1 w = — ^ WiGA . 

c(c-l) t-r 1 J J 



c(c-l) , 



Moreover 



Thus, putting together equations d!69l >. (I172I) . d!77t and (11781) yields 



with 

c-l 



(177) 



(178) 



(179) 



detF 2( JV - 1 )/ 2 2y^Fx c/2 

X = (E^XiO-i) " ■ ( 18 °) 
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Standard & Poor's 500 
10 1 , , , 




Raw Returns 



Figure 1: Graph of the Gaussianized Standard & Poor's 500 index returns versus its raw returns, during the 
time interval from January 03, 1995 to December 29, 2000 for the negative tail of the distribution. 
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a 



Figure 2: Logarithm In of the ratio of the recurrence time T to a reference time To for the recurrence 

of a given loss VaR as a function of j3 defined by (3 = VaR* (resp. VaR) is the Value-at-Risk over a 

time interval Tq (resp. T). 
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